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Using a set of hydrodynamical simulations of 9 galaxy clusters with masses in the range 
1.5 X 1O'"^M0 < Mvir < 3.4 X lO'^M,;;, we have studied the density, temperature and X-ray 
surface brightness profiles of the intra-cluster medium in the regions around the virial radius. 
We have analyzed the profiles in the radial range well above the cluster core, the physics 
[ of which are still unclear and matter of tension between simulated and observed properties, 

and up to the virial radius and beyond, where present observations are unable to provide 
any constraints. We have modeled the radial profiles between 0.3/?200 and 3/?200 with power 
laws with one index, two indexes and a rolling index. The simulated temperature and [0.5-2] 
Q ■ keV surface brightness profiles well reproduce the observed behaviours outside the core. The 

i-{ I shape of all these profiles in the radial range considered depends mainly on the activity of the 

. gravitational collapse, with no significant difference among models including extra-physics. 

■ The profiles steepen in the outskirts, with the slope of the power-law fit that changes from 

-2.5 to -3.4 in the gas density, from -0.5 to -1.8 in the gas temperature, and from -3.5 to -5.0 
. in the X-ray soft surface brightness. We predict that the gas density, temperature and [0.5-2] 

^\ ' keV surface brightness values at /?200 are, on average, 0.05, 0.60, 0.008 times the measured 

H I values at 0.37?200- At 27?200, these values decrease by an order of magnitude in the gas density 

and surface brightness, by a factor of 2 in the temperature, putting stringent limits on the 
detectable properties of the intracluster-medium (ICM) in the virial regions. 

Key words: cosmology; miscellaneous - methods: numerical - galaxies: cluster: general - 
X-ray: galaxies. 



1 INTRODUCTION 

Galaxy clusters form in correspondence of the highest peaks of the 
fluctuations in the cosmic dark matter density field and collapse un- 
der the action of gravitational attraction to define structures with a 
mean density enhanced with respect to the cosmic value by a fac- 
tor of few hundreds. They reach the virial equilibrium over a vol- 
ume with a typical virial radius that indicates, as a first approxima- 
tion, the regions where the pristine gas accretes on the dark matter 
halo through gravitational collapse and is heated up to millions de- 
grees through adiabatic compression and shocks. This intracluster- 
medium (ICM), that represe nts about 80 (15) p er cent of the bary- 
onic (total) cluster mass (e.g. lEttori et alj2004l) . becomes, then, X- 
ray emitter mainly through bremsstrahlung processes, thus allow- 
ing one to trace and characterize the distribution of the baryonic 
and dark matter components. 



From an observational point of view, the measurement of the 
properties of the ICM is possible only where the X-ray emis- 
sion can be well resolved against the background (both instru- 
mental and cosmic). At the present, data permi t to firmly charac- 



terize observables lik e surface brightness (see iMohr et aL 19991: 



lEtto ri&FabiM]|l999l) and gas temper ature (see [ Markevitch et al.l 
1998 ; iDe Griuidi & M olendi 2002; jPointecouteau et alj l200a 
Vikh linin et alj|2006l : lzhang et al. 2006) out to a fraction (~ 0.5 - 
0.6) of the virial radius. Only few examples of nearby X-ray bright 
clusters with surface brightness estimated out to the virial radius 
are available, thanks to the good spatial resolution, the large field- 
of-view an d the low instrumental background o f the ROSAT/PSPC 
instrument ( IVikhIininetal.ll 19991 : In eumann 2005). However, since 
the regions not-yet resolved occupy almost 80 per cent of the total 
volume, they retain most of the information on the processes that 
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Figure 1. From top to bottom: maps of the projected gas density, mass-weighted temperature and soft (0.5-2 keV) X-ray emission of the gl cluster (c.v/ model) 
for the 100 per cent (left) and 95 per cent (right) volume-selection scheme. The circles indicate the virial radius. The size of the side of each map is 12 Mpc, 
so they cover roughly up to 2.5i?200- 
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Table 1. List of main pliysical parameters of tlie 9 clusters taken from the ovisc runs. The objects are divided into 2 subsamples according to their virial mass. 
Note that Tvi, is the average mass-weighted temperature of the gas up to Ryir while 7200 is the temperature of the gas at ^200- 





Sample A: Afyir > lO' 






Sample B: Myi, 


< lO'^Mg 




Cluster name 


gl 


g8 




g72 


g676 


g914 


gl542 


g3344 


g62I2 


Mv„ (lO'^M.s) 


2.12 


3.39 


1.90 


1.96 


0.15 


0.17 


0.15 


0.16 


0.16 


«vir (Mpc) 


3.37 


3.94 


3.25 


3.28 


1.40 


1.46 


1.40 


1.43 


1.43 


7;,r (keV) 


6.49 


9.89 


5.78 


5.44 


1.25 


1.28 


1.18 


1.26 


1.22 


R200 (Mpc) 


2.48 


2.91 


2.39 


2.40 


1.04 


1.08 


1.04 


1.05 


1.05 


T200 (keV) 


4.18 


7.52 


4.18 


4.09 


0.82 


1.03 


0.83 


0.88 


0.84 



characterize the accretion and evolution within the cluster of the 
main baryonic component l lMolendill2004h . 

In the present work, we study the distribution of the ICM 
density, temperature and surface brightness in a set of simulated 
galaxy clusters by looking at their radial profiles in the range be- 
tween 0.3 — 3 X /?2{)(C1 We avoid considering the regions where the 
clusters evolve a core with active cooling and suspected feedback 
regulation, being this balance still matter of tension between the 
present si mulated constraints and the observed counterparts (e.g. 
iBorgani et a l. 2004). Therefore, we concentrate our analysis on 
cluster volumes where the gas density and temperature suggest that 
radiative processes are not dominant and the main physical driver 
is just the gravitational collapse, that is expected to be well de- 
scribed in the present numerical N-body simulations. We consider 
objects extracted from a dark matter only cosmological simulation 
that were resimulated at much higher resolution by including the 
gas subjected to gravitational heating and other physical treatments 
like cooling, star formation, feedback, thermal conduction and an 
alternative implementation of artificial viscosity. 

We organize this paper in the following way: in the next Sec- 
tion, we describe our dataset of simulated clusters and their differ- 
ent physical properties; in Section[3]we discuss the results on the 
slopes of the gas temperature, density and surface brightness pro- 
files in the outskirts (i.e. r > 0.3/?20o) of the simulated objects and 
compare our results with the present constraints obtained from X- 
ray observations of bright clusters. We summarize and discuss our 
results in Section|4l 



2 THE SIMULATED CLUSTERS 

For the purpose of this work we use a set of 9 resimulated galaxy 
clusters extracted from a dark matter (DM)-only simulation of a 
"concordance" ACDM model with flm = 0.3 and flyv = 0.7, using 
a Hubble parameter h = 0.7 (the Hubble constant being Hq = h 100 
km s^'Mpc^^) and a baryon fraction of Q.\, = 0.0l9h^^. The ini- 
tial conditions of the parent simulation, that followed the evolution 



In this paper we define the typical physical quantities of clusters as fol- 
lows: i?2oo (R500) is the radius of a sphere centered in a local minimum of 
the potential and enclosing an average density of p = 200(500)pci , with Pci 
being the critical density of the Universe. The virial quantities are defined 
using the density threshold obtained by the spherical collapse model with 
the cosmological parameters of the simulations. So the virial rad ius /fyj, 
encloses an average density of pv„ ^ 104pcr (see lEke et alll 19981 eq. 4). 
Consequently the virial mass Mvip is the total mass (DM plus baryons) of 
the cluster up to Ryji . 



of a comoving box of 419h^ Mpc per size, were set consider- 
ing a cold da rk matter power spe ctrum normalized by assuming 
ag = 0.9 (see lYoshida et al.ll200lh. O ur simulations were carried 
out using GADGET-2 (Springel 2005), a new version of the Tree- 
Smoothing Par ticles Hydrodynamics (SPH) parallel code GADGET 
dSpringelet^ boOl ) which in cludes an entr opy-conserving for- 
mulation of SPH cSpringel & Hernquislll2002h and a treatment of 
many physical processes affecting the baryon component that can 
be turned on and off to study their influence on the thermodynamics 
of the gas. In the volume of the parent simulation 9 different haloes 
were identified, spanning different mass ranges: 4 of them with the 
typical dimension of clusters and 5 of small clusters or group-like 
objects . Then using the "Z oomed Initial Conditions" (ZIC) tech- 
nique ( lTormenetal] |l997) they were resimulated by identifying 
their corresponding Lagrangian regions in the initial domain and 
populating them with more particles (of both DM and gas), while 
appropriately adding high-frequency modes. At the same time the 
volume outside the region of interest was resimulated using low- 
resolution (LR) particles in order to follow the tidal effects of the 
cosmological environment. The setup of initial conditions of all 
the resimulations was optimized to guarantee a volume around the 
clusters of --^ 5i?vir free of contamination from LR particles. This 
was obtained using an iterative process as follows. Starting from a 
first guess of the high-resolution (HR) region that we want to res- 
imulate, we run a DM-only re-simulation. Analyzing its final out- 
put, we identify all the particles that are at distances smaller than 5 
^vir from the cluster centre in order to identify the corresponding 
Lagrangian region in the initial domain. Applying ZIC, we gen- 
erate new initial conditions at higher resolution and run one more 
DM-only resimulation. This procedure is iteratively repeated until 
we find that none of the LR particles enters the HR region, which 
could be possible because of the introduction of low-scale modes. 
Thus we can safely say that these resimulations are representing 
the external regions of our cluster without any spurious numerical 
effect. 

These resimulations are set to have a mass resolution of 
'"DM = L13 X 10^/?" 'M0 for DM particles and mgas = 1.69 x 
IO'^H^^Mq for baryons, so that every resimulated cluster is re- 
solved with between 2 x 10^ and 4 x 10^ particles (both DM 
and gas) depending on its final mass. All the simulations have 
a (Plummer-equivalent) softening-length kept fixed at e = 30h^^ 
kpc comoving at z > 5 and switched to a physical softening length 
£ = 5/7^' kpc at lower redshifts. 

In order to study the impact of different physical processes on 
the clusters properties, we simulated every cluster using 4 different 
sets of physical models: 

(i) Gravitational heating only with an ordinary treatment of gas 
viscosity (which will be referred as the ovisc model). 
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(ii) Gravitational lieating only but with an alternative implemen- 
tation of ar tificial viscosity {Ivisc mod el), following the scheme 
proposed bv Morris & MonaghanI l ll997h . in which every gas parti- 
cle evolves with its own viscosity parameter. With this implemen- 
tation the shocks produced by gas accretion are as well captured as 
in the ovisc model, while regions with no shocks are not character- 
ized by a la r ge res idual numerical viscosity. As already studied in 
iDolag et al.l |2005^, this scheme allows to better resolve the turbu- 
lence driven by fluid instabilities, thus allowing clusters to build up 
a sufficient level of turbulence-powered instabilities along the sur- 
faces of the large-scale velocity structures present in cosmic struc- 
ture formation. 

(iii) Runs which implement a treatment of cooling, star for- 
mation and feedback (csf model). The star formation is followed 
by adopting a sub-resolution multiphase model for the interstellar 
medium which includ es also the fee dback from supemovae (SN) 
and galactic outflows jSpringel & H ernquist 2003). The efficiency 
of SN to power galactic winds has been set to 50 per cent, which 
turns into a wind speed of 340 km/s. 

(iv) Runs like csfhut also including th ermal conduction jcsf c 
model), adopting the scheme described by ljubelgas et alj ( |2004) . 
This implementation in SPH, which has been proved to be stable 
and to conserve the thermal energy even when individual and adap- 
tive time-steps are used, assumes an isotropic effective conductivity 
parametrized as a fixed fraction of the Spitzer rate ( 1/3 in our simu- 
lations). It also accounts for saturation which can become relevant 
in low-density regions. For more detail s on the properties of simu- 
lated clusters using this model, see also lDolag et alj ( [2004) . 

The masses, radii and temperature of the 9 clusters for the 
ovisc model are summarized in Table [T] The other models have 
similar results for masses and radii with variation of few percent 
between the models while temperatures have larger variations with 
csf and csfc models having slightly higher temperatures (~ 8 per 
cent). In our following analyses we will separate the clusters into 
two different samples: sample A (4 clusters) and B (5) according to 
their virial mass being larger or smaller than IO'^Mq, respectively. 

When considering the environmental properties of these 
haloes, we must point out that the 4 objects of sample A are among 
the most massive clusters of the volume of the parent simulation 
while the 5 groups were chosen at random with the only condition 
to be far away from the 4 clusters. At the end of the resimulations 
all the haloes of sample B appear to be isolated, e.g. there is no 
halo with mass M > IO'-^'^A^'Mq at distance lower than 5i?vii-; on 
the contrary the big clusters have several small structures inside the 
volume of their resimulation and all of them underwent significant 
merging activities. As we will see in Section [S] this aspect is im- 
portant to understand the different properties of the X-ray profiles 
of the two samples. 



^ Notice that in this paper we prefer to quote all the temperatures by us- 
ing the mass-weighted estimator because it is more related to the energetics 
involved in the process of stru c ture formation. As s hown in earlier p apers 
(see, e.g. iMazzotta et al.ll2004 iGaidini et ai]|2004 liasia et al.ll200a) , the 
application of the emission-weighted temperature, even if it was originally 
introduced to extract from hydrodynamical simulations values directly com- 
parable to the observational spectroscopic measurements, introduces sys- 
tematic biases when the structures are theiTnally co mplex. We verifie d that 
the use of the spectroscopic-like temperature of Ma zzotta et alj J2004) does 
not produce significant changes in our results for kgT >\ keV. 



2.1 Preparation of the dataset 

Focusing our attention on the cluster outskirts, we do not consider 
in our study the effects of the cluster core on the overall proper- 
ties, avoiding to deal with the theoretical, numerical and observa- 
tional uncertainties on the physics of the X-ray emitting plasma 
in the cluster central regions. Furthermore we must take into ac- 
count that the mass enclosed in a given radial bin in the external re- 
gions of clusters (r > -R200) is dominated by the presence of dense 
subclumps of the typical dimension of galaxies or small galaxy- 
groups (M < 1O'"'M0). Evidences of the existence of these clumps 
in the regions inside i?200 have also been shown by X-ray observa- 
tions, but it is still not clear whether the frequency and the emission 
properties of the simulated clumps agree with the observed ones: 
this is due both to the low surface brightness of these objects that 
makes them difficult to resolve and to uncertainties in the theoret- 
ical modeling (e.g. simulations including cooling tend to produce 
more clumps and simulations with feedback from SNe end up with 
smaller clumps than non-radiative ones). 

Anyway, since (i) we are interested only in the average prop- 
erties of galaxy clusters and (ii) these small clumps are generally 
masked out in observational works, a direct comparison of our 
simulations with theoretical and observational results requires the 
introduction of a method to exclude these cold and dense clumps 
from our analysis. It is quite difficult to define a precise and unique 
criterium to identify these condensed regions in the simulation 
volume from their dynamical and thermodynamic al properties, 
therefore we decided to identify the volumes corresponding to 
non-collapsed regions. For this purpose we proceed in the follow- 
ing way: for every shell that corresponds to a radial bin of the 
profile, we sort the gas particles in decreasing volume order (i.e. 
starting from the more diffuse particles and going up to the denser 
ones, being the volume associated to the /-th particle defined as 
Vi = OTgas/Pgas.i) then we compute the profiles calculating the 
mean of a given quantity considering only the particles summing 
up to a given percentage of the total volume of the shell. In order 
to identify the physical properties of the regions excluded by our 
analysis, we plot in Fig.[T]as an example the maps of the projected 
density, mass-weighted temperature and soft (0.5-2 keV) X-ray 
surface brightness of the gl cluster considering (i) all the particles 
and (ii) those that are within the 95 per cent volume of each radial 
shell. We can clearly see from these maps that the volume- selective 
method actually separates the cold and dense clumps from the rest 
of the object. Anyway these maps show that this criterion is much 
more effective in excluding dense regions than cold ones. However, 
since the X-ray emission is strongly dependent on the density of 
the gas, we can safely say that it well mimics the observational 
technique of masking the bright isolated regions. 

We also study the shapes of X-ray surface brightness profiles 
of our simulated clusters. Since we are interested in the outskirts 
of galaxy clusters we cannot adopt a pure bremsstrahlung emis- 
sion model because in the corresponding temperature ranges (T < 
I keV) line emission gives a significant contribution. Therefore in 
order to calculate the profiles we adopt a MEKAL emission model 
jMewe et alJll985l : lLiedahI et alj[T99a) assum ing a gas metallicity 
Z = 0.2Zq (using the solar value tabulated in ' Anders & Grevessd 
constant for every particle of the simulation: although this 
value c an be an underestimate o f the metallicity of the clusters' cen- 
tre (see lDe Grandi et ^|2004| . for more detail), the chosen model 
can be considered a good representation of the spectral properties 
of the clusters' external regions. In the following analyses we con- 
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Figure 2. Radial profiles for different quantities as estimated for the gl cluster. From top to bottom: gas density (normalized to the cosmic mean density, 
< p >= ilbPc), mass-weighted temperature and X-ray surface brightness for the soft (0.5-2 keV) and the hard (2-10 keV) band calculated considering the 
clusters at z = 0. Left column: ovisc model. Right column: ci/' model. These profiles are obtained from different prescriptions: the black line shows the average 
of all the gas particles, the orange lines are ave rages of the different volume percentages from 10 to 90 per cent and the green lines are from 91 to 99 per cent. 
The extragalactic unresolved background from lHickox & Markevitchl j2006l) in the [0.5-2] keV and [2-8] keV band of is indicated by the shaded region in the 
soft and hard X-ray surface brightness panels, respectively. 
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sider the objects at redshift z = 0. To compute the surface bright- 
ness profiles, we adopted the same volume-selection scheme of the 
density and temperature profiles with the only difference that, since 
we want to obtain the surface brightness profiles of our objects, we 
need to take into account the two-dimensional projection on the sky 
plane. To this purpose, we consider circular ring sections instead of 
spherical shells: for every ring, we select all the particles that fill 
a given volume fraction /, we calculate their luminosity and we 
obtain L j by summing all of them. Finally we normalize it to take 
into account also for the remaining volume and use the normalized 
luminosity L,-,„g = y- to calculate the surface brightness (see the 
effects of this method on the [0.5-2] keV X-ray surface brightness 
shown in the bottom panel of Fig.[T](. This method works efficiently 
for subtracting almost all the relevant substructures in the surface 
brightness profiles, except for two massive subhaloes present in the 
gl cluster simulation at a distance of 2-3 /?2()0 from the center; only 
for this cluster we excluded two separate angular regions (45 de- 
grees wide each) from the computation of its X-ray surface bright- 
ness profiles. 

In Fig. |2l we show the gas density, temperature and surface 
brightness profiles into two different bands for a massive clus- 
ter (gl) assuming two different physical models: the "ordinary- 
viscosity non-radiative" (ovisc) and the "cooling + star formation + 
feedback" (csf) one. We plot the average profiles together with the 
ones obtained with different cuts in volume. The profiles in Fig.|2] 
show that eliminating the fraction of mass concentrated in 1 per 
cent of the volume is sufficient to make the gas density, temperature 
and surface brightness profiles much more regular than the average 
ones that consider all the particles in the shell. For these reasons 
in the following analyses, we will mainly concentrate in fitting the 
99 per cent-volume profiles and discuss the results obtained by cut- 
ting the volume at different percentages in Appendix lAl By the way 
we also note that the temperature and density profiles obtained with 
this method for the ovisc model well reproduce the shape of the the- 
oretical profiles obtained by Ascasibar et al. (2006), if we assume 
a polytropic index y ~ 1 . 1 8 — 1 .20. 

In Fig. [3] we compare the 99 per cent volume profiles for one 
of the largest and the smallest cluster of the sam ples using the 
four d ifferent physical models. As already noted in Bor gani et al.l 
( l2004h . radiative cooling and star formation selectively remove 
low-entropy gas from high-density regions. Consequently mod- 
els with cooling tend to produce clusters less dense at the centre 
(r < 0.1i?20o)' in these regions we find that the gas density in non- 
radiative models is higher by a factor of ~ 3.5 in the most massive 
clusters and more than 5 in less massive ones. The same cooling 
processes generate a lack of pressure near the centre of the clusters 
that heats the infalling gas up to a factor of ~ 1 . 8 hi gher than in non- 
radiat ive simulations for all mass scales (see also iTornatore et al] 
l200j ). As a result, the X-ray emission in the internal regions is sig- 
nificantly reduced in the csf and csfc models. The external regions 
of the clusters are less affected by the different physics because the 
low-density gas is less prone to cooling and feedback effects. As 
we will show later, this makes the results on the external regions of 
clusters only slightly d ependent on the phy sical model adopted in 
the simulation (see also lRomeo et al ] |2006l , for influences of other 
physical processes). 

In Fig.|4]we compare the 99 per cent volume profiles for the 
4 most massive clusters of the dataset (sample A). The density pro- 
files have a very similar shape up to about 2^200 indicating a uni- 
versal function. The regularity of the 99 per cent volume profiles 
extends also for distances > 3/?200- The temperature profiles are 
more irregular (note that the temperatures are normalized using 



T200, the temperature at R200 )'■ the bumps present in the gl pro- 
file near the centre and in g72 at r ~ 1 .7/?200 are an effect of shocks 
at r ~ 0.6/?200 ™d ~ 2.1^200> respectively, due to recent major 
mergers. 

In Table[2]we show the ratios between the values of the 99 per 
cent volume profiles at different radial distances and their values 
at 0.3i?vii- obtained with the ovisc model. The regularity of these 
profiles can be seen by the low dispersion between the values of 
the clusters of the same sample. 



3 RESULTS ON THE OUTER SLOPES OF THE RADIAL 
PROFILES 

To describe the behaviour of the radial density and temperature pro- 
files in the cluster outskirts, we adopt a broken power-law relation 
given by the expression 



.-hi 



ifx^Rh/R2oo 
if X > Rb/Rim 



(1) 



in its logarithmic form, where r is the distance from the centre, 
X = r//?200 ™d a, b\, b2 and Rj, are the free parameters, repre- 
senting the normalization, the inner slope, the outer slope and the 
radius at which the slope changes, respectively. The fit is computed 
in the interval 0.3 ^ x ^ 2.7. We fit the profiles for the 9 clusters 
separately and then calculate the average and dispersion for every 
parameter. The best fit relations for the ovisc model are shown to- 
gether with the profiles of single clusters in Fig.|5] 

The gas density profiles show a change in the slope of the 
power law at about 1.1-1.3 i?2(K)i with an internal slope of 2.4-2.5 
with no significant dependence on the physical model or the volume 
selection scheme adopted (see the Appendix for details on the effect 
of different volume cuts) and also a very low dispersion between 
the different objects. This indicates that the behaviour of the density 
profiles is well determined from outside the core to more than i?200- 
In the outskirts the slope increases up to b2'~^3.4 for all the models 
showing more dispersion between the 9 haloes due to the fact that 
in these regions the influence of the environment is becoming more 
important with respect to the potential well of the cluster itself. 

In the temperature profiles, the steepening towards the exter- 
nal regions is more prominent and the cutoff is around r = 1.3 — 
1.5i?200 with significant dispersion between the clusters. The best 
fit internal slope is significantly affected by the physics: the profiles 
obtained with non-radiative simulations have slopes of hi ~ 0.4 
while csf and csfc models produce steeper profiles with slopes of 
bi ~ 0.6. Also in the outskirts the non-radiative models tend to 
produce slightly shallower profiles with the difference less signifi- 
cant due to the higher dispersion between the different haloes. The 
different slopes of the profiles is an effect of the cooling that creates 
an higher gradient of temperature as can be seen in Fig. [3] 

In Fig. |2] we also show the soft (0.5-2 keV) and hard (2-10 
keV) X-ray surface brightness profiles of the gl object for the dif- 
ferent physical models. These plots indicate that models with cool- 
ing produce a higher number of clumps that dominate the soft X- 
ray emission in the external regions. These clumps are indeed those 
that are well eliminated by our selection in volume. We also obtain 
that the cluster surface brightness at z = r eaches the value of the 
unresolved X-ray background estimated by iHickox & Markevitchl 
1 I2OO6I) of 1.76±0.32 in the [0.5-2] keV band and 3.5± 1.7 in the 
[2-8] keV band in units of 10^'^ erg s^' cm^^ deg^^, at approxi- 
mately 0.8i?200 and OAR200 in the soft and hard band, respectively; 
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Figure 3. From top to bottom: comparison between the gas density (normalized to the cosmic mean density, < p >= iltpc), mass-weighted temperature, soft 
and hard X-ray surface brightness profiles extracted from the 99 per cent vo lume of the clusters (left column) and g6212 (right column) simulated by using 4 
different physical models. A dashed line indicates the functional form from lvikhlinin et al] j2006l eq. 9) that well reproduces the behaviour of the temperature 
profile of nearby bright galaxy clusters observed with Chandra. The functional form is normalized to the average values of the temperature obtained from the 
4 models at 0.5i?200- The extragalactic unresolved background from lHickox & MarkevitclJ j20O6l) in the [0.5-2] keV and [2-8] keV band of is indicated by 
the shaded region in the soft and hard X-ray surface brightness panels, respectively. 
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Figure 4. Comparison between the 99 per cent volume profiles of density (normalized to the cosmic mean density, < p >= ilbpc) (top panels) and mass- 
weighted temperature normalized to Tim (lower panels) for the clusters gl (red), g8 (blue), g51 (green), g72 (orange). Left and right columns refer to the ovisc 
and cs/ model, respectively. 



Table 2. Ratios between the profile values at different radial distances and their values at 0.3R2oo- We show the average and the standard deviation between 
the two cluster samples for the ovisc model. 



Radial distance (R200) 


0.3 


0.5 


0.7 




1.0 




2.0 






3.0 
















Sample A 
















P 


1.0 


0.31 ±0.02 


0.127 ±0.004 




(4.8 ±0.2) X 10" 


-2 


(5.6±0.8) > 


( 10" 


-3 


(1.7±0.5) > 


< 10" 


-3 


T 


1.0 


0.88±0.06 


0.72 ±0.06 




0.60 ±0.06 




0.34 ±0 


09 




0.16±0 


03 




■^S|o,5_2|kcV 


1.0 


0.17 ±0.03 


(4.0±1.1) X 10- 


-2 


(8.4±1.2) X 10" 


-3 


(2.3±0.9) > 


< 10" 


-4 


(6.3±4.3) > 


< 10" 


-5 


■5'S|2-10|kcV 


1.0 


0.15 ±0.02 


(3.2±0.1) X 10" 


-2 


(5.5±0.1) X 10" 


-3 


(7.1±1.8) > 


( 10" 


-5 


(6.4±5.7) > 


< 10" 


-6 












Sample B 
















P 


1.0 


0.31 ±0.02 


0.13±0.02 




(4.8 ±0.7) X 10" 


-2 


(6.0±0.2) X 10" 


-3 


(1.6±0.2) > 


< 10" 


-3 


T 


1.0 


0.84 ±0.04 


0.72 ±0.04 




0.58 ±0.03 




0.29 ±0 


06 




0.14±0 


02 




'^^[0.5-2|keV 


1.0 


0.16 ±0.02 


(4.2±1.1) X 10" 


-2 


(7.9 ±1.0) X 10" 


-3 


(1.8±0.9) > 


c 10" 


-4 


(1.5±1.0) > 


< 10" 


-5 


SB|2-10|keV 


1.0 


0.11 ±0.01 


(2.1 ±0.7) X 10" 


-2 


(2.3 ±0.4) X 10" 


-3 


(1.7±1.4) > 


( 10" 


-5 


(5.6±7.2) > 


< 10" 


-7 



considering the effe ct of redshift dimmin g and that, according to 
the recent results of iRoncarelli et al] feOOffl . a significant fraction 
of this background can be due to the diffuse gas of non-virialized 
objects, the perspective of observing the regions around the clusters 
virial radius is extremely challenging. 

Since we observe a steepening of the slope of the X-ray pro- 
files going to the external regions, we fit them with a rolling-index 
power law function 



where x = r/R2oo, and a, iin,;,, and bj^ax are the free parameters, 
with fc^in and bj^^x representing the slope in the two asymptotic 



cases, x~fO and x S> 1, respectively. We fit the two samples sepa- 
rately in the interval 0.3 ^ r/R2oo ^ 2.7. The results for the 4 phys- 
ical models are quoted in Tab. |3] The average best-fit values on the 
slopes of the surface brightness of the objects simulated with the 
ovisc model are shown in Fig.|6] Overall, the fits well reproduce the 
shape of the surface brightness in the interval 0.3 ^5 y/RiOQ ^ 2.7. 
Our sample of galaxy "groups" (sample B) presents a soft X-ray 
profile that is (i) shallower in the centre (^rnin.B ~ 2.7, iimin.A ~ 3.4) 
and (ii) steeper in the outskirts (fcmax.B ~ 7.3, fcmax.A ~ 5.5) than 
the one estimated in the sample of massive clusters (sample A). 
The following reasons can be considered to explain this behaviour: 
(i) the ICM temperature of the objects in sample B at R2Q0 is be- 
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Figure 5. Profiles (99 per cent volume) of the 9 clusters obtained with the ovisc model. The orange lines are for the 4 clusters of sample A, the green lines are 
for the 5 groups of sample B. We fit the data with a broken power-law relation of eq.[T]in the interval 0.3 ^ r/i?200 ^ 2.7 and plot (black line) the relation 
assuming the average of the best-fit parameters between the 9 objects. 




Figure 6. X-ray surface brightness profiles (99 per cent volume, [0.5-2] and [2-10] keV band) of the 9 clusters for the ovisc model. The orange lines are 
for the 4 clusters of sample A, the green lines are for the 5 groups of sample B. We fit the data of the two samples separately with a rolling-index power-law 
relation of eq.[2]in the interval 0.3 ^ r/R2oo ^ 2.7 and plot (black lines) the relations assuming the average of the best-fit parameters between the objects of 
the samples. 



low 1 keV (see Table [TJ and decreases at r > i?200 to values that 
becomes comparable to (and less than) the lower end of the [0.5-2] 
keV band here considered, resulting in a cut-off of the soft X-ray 
signal. This also explains the steeper slope of the hard X-ray pro- 
files in the less massive objects; (ii) the higher cooling efficiency 
in groups, with respect to rich clusters, enhances the selective re- 
moval of low-entropy gas from the hot phase, thus producing shal- 
lower profiles of X-ray emitting gas; (iii) as discussed at the end 
of Section|2l while our clusters are the most massive objects in the 
volume of the parent simulation and are still accreting material, the 
"groups" appear to be isolated, thus with less evidence of recent ac- 
cretion activity. This results in slightly shallower profiles near the 
centre. 



3.1 Comparison with observational constraints 

Our results on the external shape of the gas density, temperature 
and surface brightness profiles in simulated galaxy clusters can be 
compar ed with recent es timates obtained for nearby X-ray bright 
objects. i Neumani] ( l2005h discusses the outer slope of the soft X- 
ray surface brightness in 14 clusters observed with ROSAT/PSPC. 
In order to compare our results with the ones of that work, we fit the 



99 per cent volume X-ray profiles in different intervals for the two 
clusters samples with a single power-law function, S{x) = ax^'', 
where a and b are the two free parameters and x = r/i?200- We 
avoid to consider the interval 0.1 ^ x ^ 0.3, owing to the fact that 
this region is very model-dependent and cannot be represented by a 
single power-law relation. Our results are listed in T able [4] together 
with the corresponding results (iobs) of iNeumannI |2005). In the 
interval 0.3 ^ x ^ 1.2, for the sample of clusters labeled "4" that 
better matches the virial temperatures of our sample A. Neumannl 
( I2OO5I) measures a slope of 3.79 ± 0.38 in good agreement with our 
average values of (3.78,4.03), being the lower values measured in 
objects simulated with the inclusion of extra-physics (models csf 
and csfc). In the interval 0.7 ^ x ^ 1.2, the observational constraints 
are looser and are between 5.73^} for the 7 most massive clus- 
ters in the Neumann's sample and 7.22['[} gg for the whole dataset. 
Our constraint of ~ 4.3 — 4.5 lays on the lower end of this distribu- 
tion, but still just lo away from the results obtained for the hottest 
sample. The sample B shows profiles with mean slopes consistent 
with those obtained in sample A within the measured dispersion 
(see Table|4l(. 

For what concerns the temperature profile, a comparison can 
be done with the functional form that reproduces the behaviour 
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Table 3. Distribution of tlie best-fit values of the rolling-index power-law functions (see eq.|2) of the soft (0.5-2 keV) and hard (2-10 keV) X-ray surface 
brightness profiles separated between the two samples A and B and for different physical models. The profiles are extracted from the 99 per cent volume- 
selection scheme and fitted over the radial range 0.3 ^ r/Rim ^ 2.7. 



Soft X-ray band 
[0.5-2] keV 
Sample A Sample B 





h 

''mill 






^max 


ovisc 


3.45 ±0.39 


5.67 ±0.67 


2.82±0.38 


7.21 ±1.43 


Ivisc 


3.36±0.33 


5.73 ±0.80 


2.73 ±0.79 


7.72 ±1.93 


csf 


3.49±0.21 


5.33 ±0.57 


2.50±0.28 


7.14±1.14 


csfc 


2.76 ±0.28 


6.36±1.09 


1.99 ±0.50 


7.80 ±1.76 






Hard X-ray band 








[2- 


-10] keV 






Sample A 


Sample B 




^mtii 






^inax 


ovisc 


3.10±0.55 


7.97 ±1.01 


3.42±0.71 


9.82 ±2.50 


Ivisc 


3.08 ±0.29 


7.87 ±0.90 


3.36±1.72 


10.48 ±5.26 


csf 


3.25 ±0.41 


7.56 ±0.73 


3.25 ±0.43 


9.69 ±2.66 


csfc 


1.94 ±0.36 


9.31±2.11 


2.87±1.71 


10.95 ±5.79 



Table 4. Distribution of the best-fit values of the power-law slopes of the X-ray surface brightness profiles separated between the two samples A and B (b^ and 
hs, respectively), the two X-ray bands and different physical models. The profiles extracted from the 99 per cent volume-selection scheme and fitted over the 
different radial ranges A',^i„ —x^^y, (in units of R200-) Values of io^s the observational constraints in the same radial range from lNeumannI i2005h (see text 
in Section im . 



Soft X-ray band Hard X-ray band 

[0.5-2] keV [2-10] keV 





-^'min 


-^max 


lu 


be 


bobs 


bA 


be 


ovisc 


0.3 


- 1.2 


4.05±0.03 


4.09±0.01 


3.79±0.38 


4.46±0.05 


5.23±0.03 


Ivisc 


0.3 


- 1.2 


4.04±0.04 


4.14±0.09 




4.47±0.05 


5.32±0.17 


csf 


0.3 


- 1.2 


3.99±0.00 


3.81 ±0.02 




4.50±0.02 


5.07±0.04 


csfc 


0.3 


- 1.2 


3.79±0.02 


3.57±0.04 




4.19±0.03 


5.03±0.12 


ovisc 


0.7 


- 1.2 


4.49±0.19 


4.76±0.32 


-'■'-'-1.26 


5.15±0.41 


6.42±0.76 


Ivisc 


0.7 


- 1.2 


4.54±0.32 


4.66±0.19 


5.26±0.44 


6.20±0.75 


csf 


0.7 


- 1.2 


4.31±0.28 


4.45±0.11 




5.18±0.37 


5.92±0.13 


csfc 


0.7 


- 1.2 


4.29±0.17 


4.35±0.14 




4.89±0.24 


6.42±0.61 


ovisc 


1.2 


-2.7 


5.15±0.62 


6.37±0.93 




6.95±0.88 


8.43±1.70 


Ivisc 


1.2 


-2.7 


5.21±0.59 


6.77±0.95 




6.88±0.54 


9.01±2.72 


csf 


1.2 


-2.7 


4.90±0.53 


6.23±0.67 




6.70±0.74 


8.40±1.52 


csfc 


1.2 


-2.7 


5.67±0.63 


6.58±0.88 




8.14±1.32 


9.14±3.36 



of the deprojected X-ray temperature profile at r > 0.05i?2()0 of 
1 3 low-redshift clust ers observed with Chandra as presented in 
IVikhlinin et all bOOd see their equation 9): 



r(d)< 



1 



(d/0.045)''^+0.45 
(d/0.045)i-9 + l ■ "(r+ (d/0.6)2) 



2U).45 ' 



(3) 



where d = r/R^Q^). We overplot this function normalized to the av- 
erage temperature value at 0.5i?200 to our profiles in Fig.[3l 

While it is known that in general hydrodynamical simulations 
cannot reproduce the steepening of the observed temperature pro- 
files at the centre (see lBorgani et al.ll2004l I2OO6I for more detail), 
we can see that the agreement in the external regions is good for the 



mo st massive clusters, as a l ready noted bv lMarkevitch et all fT99^ 
and lDe Grandi & Jvlolendil ilOOA The non-radiative models show 
a good agreement with this form, so that when we normalize the 
function at 0.5i?200 the difference between the values quoted in Ta- 
ble|2]is less than I0 up to 2i?200 for both clusters and groups. We 
only find small deviations from the observed profile that become 
more significant at r < 0.3i?200 for simulations including cooling 
and at r > R200 where non-thermalized accreting material is domi- 
nant. The latter effect is particularly evident in low-mass systems. 
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3.2 Implications on X-ray properties of the cluster virial 
regions 

Using our profiles we can predict the steepness of the profiles 
in the external regions of galaxy clusters. We fit the soft (0.5- 
2 keV) and hard (2-10 keV) X-ray profiles also in the interval 
1.2 ^5 r/RiOQ ^ 2.7. The results are shown in Table|4] Our profiles 
are steeper in the external regions for all the models and for both 
the cluster samples: the steepening is more evident in sample B 
suggesting a break of self-similarity at the scales around 2i?20()- As 
already noted at the end of SectionjS] this is mainly due to the fact 
that the temperature in the external regions of the haloes of sample 
B drops below 0.5 keV resulting in a cutoff of the soft X-ray sig- 
nal. From the fit with a rolling power-law, the values of t'max can 
be associated to the asymptotic behaviour at r S> Riqo- We measure 
a slope of the surface brightness that ranges from about ~ 7.5 in 
grou ps to ~ 5.8 in clusters. In terms of t he |3 value of the (3— model 
(e.g. ICavaliere & Fusco-Femianolll978l) , P fi; (1 +fcmax)/6, it ap- 
proaches an estimate between 1 .4 in groups and 1 . 1 in clusters, the 
latter being a 40 per cent s teeper than the observed values at --^ i?20() 
(IVikhlininp,t"Ii]|l999l : lNeumann..2005i) . 



4 SUMMARY AND CONCLUSIONS 

Using a set of hydrodynamical simulations, performed with the 
Tree-l-SPH code GADGET-2, composed by 9 galaxy clusters cov- 
ering the mass range 1.5 X IO'^Mq < Mvir < 3.4 X IO'^Mq and 
adopting 4 different physical prescriptions with fixed metallicity, 
we have studied the behaviour in the outer regions of the gas den- 
sity, temperature, soft (0.5-2 keV) and hard (2-10 keV) X-ray sur- 
face brightness profiles paying particular attention to the logarith- 
mic slope. Outside the central regions, where a not-well defined 
heating source can partially or totally balance the radiative pro- 
cesses, the physics of the X-ray emitting intracluster plasma is ex- 
pected to be mainly driven by adiabatic compression and shocks 
that take place during the collapse of the cosmic baryons into ac- 
creting dark matter halos. These processes are properly considered 
in cosmological numerical simulations here analyzed, allowing us 
to investigate on solid basis the expected behaviour of the cluster 
outskirts. 

Our main findings can be summarized as foUowfl 

(i) the behaviour of the profiles in the external regions of clusters 
{r > i?20o) does not depend significantly on the presence or absence 
of cooling and SN feedback. Only the inclusion of a model of ther- 
mal conduction can introduce non-negligible changes in the values 
of the slopes of the profiles. This shows that the mean behaviour of 
the X-ray emitting plasma in the cluster outskirts is mainly due to 
the gravitational force. 

(ii) The gas density profile steepens in the outskirts changing the 
slope of the power-law from ~ 2.5 to ~ 3.4 at about 1.2i?200- The 
temperature profile varies the slope of the power-law like functional 
form from 0.4 in non-radiative models and 0.6 in models with cool- 
ing and feedback, to about (1.7 — 1.9) at ~ 1.3 — 1.6i?200- 

(iii) The surface brightness in the outskirts profiles appear to be 
dominated by the presence of subclumps of cold and dense material 
that we exclude from our analyses. After removing these objects, 
the surface brightness in the [0.5-2] keV presents a profile with a 

' Note that, as in the rest of the paper, the values of the slopes here quoted 
are in fact negative slopes. 



slope that steepens towards the external regions. Fitting the profiles 
with a power-law function in different radial ranges between 0.3 
and 2.7 R2Q0 we obtain slopes between 4 and 5.5 in the 4 most mas- 
sive systems and between 3.8 and 6.5 in the 5 objects with masses 

(iv) We introduce a new rolling-index power-law functional (eq. 
O which fits well the X-ray surface profiles in the interval 0.3 ^ 
'■/^200 ^ 2.7. From the results of these fits, we obtain that our sam- 
ple of groups of galaxies (sample B) has soft X-ray profiles shal- 
lower in the centre (fcmin.B ~ 2.6, imin.A 3.3) and steeper in the 
outskirts (fcmax.B ~ 7.4, fcmax.A ~ 5.8) than our sample of clusters 
(sample A). This is mainly an effect of the drop out from the pass 
band of the soft X-ray emission, due to the low ICM temperature 
of the objects in sample B at the virial radius. 

(v) We compare our results on the shape of the temperature and 
soft X-ray surface brightness profiles with the present observational 
constraints. By using the functional form of Vikhlinin et al. (200^ 
(eq. 9) extrapolated beyond R^oo^ we verify that there is a good 
agreement with the simulated profiles over the radial range 0.3 — 
2i?200- On the surface brightness profile in the [0.5-2] keV band, 
we measure a slope in the outskirts that is consistent within la with 
the values estimated from ROSAT/PSPC data in Neumann (2005). 

(vi) From the overall shape of the simulated profiles, we pre- 
dict that (see Table O (1) the ICM temperature profile decreases 
with radius reaching at i?200 average values of 0.60 times the value 
measured at 0.3 x i?200> (2) the estimated gas density at ^^200 is 
4.8 X lO^-' times the gas density measured at 0.3i?200' (3) the sur- 
face brightness in the soft band is about 10~'^ erg s~' cm~^ deg^'^ 
at i?200 at z=0, that is a factor of few below the e xtragalactic unre- 
solved X-ray background (see lHickox & Markev itch 2006). 

This last item shows that the perspective to resolve the X-ray 
emission around the cluster virial regions is extremely challenging 
from an observational point of view, even at low redshift, unless 
tight constraints on all t he components of the measured background 
can be reached (see also lMolendil2004l) . Considering the similarity 
of the surface brightness profiles of the objects of the same mass, 
the observational technique of stacking the images of different clus- 
ters (e.g. Neumann 2005) can help in enhancing the statistics of the 
profiles in the outer regions. 

Overall, our results show how the study of simulated objects 
can give indications on the requested capabilities of future exper- 
iments in order to resolve the cluster virial regions. On the other 
hand, since the behaviour of the slopes have shown to be mainly in- 
dependent of the physical models adopted in our simulations, any 
discrepancy between our results and observations of the clusters 
virial regions could provide useful hints on the presence of addi- 
tional physics (e.g. turbulence, magnetic fields, cosmic rays) influ- 
encing the thermodynamics of the ICM. 

Finally, the work presented here further demonstrates the com- 
plementarity between numerical and observational cosmology in 
describing the properties of highly non-linear and over dense struc- 
tures. 
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APPENDIX A: EFFECTS OF VOLUME CUTTING IN 
DENSITY AND TEMPERATURE PROFILES 

In this appendix we will briefly discuss the effects of the volume- 
selection scheme adopted in this paper. We show in Table IaTI and 
IA2I the results for the fits of the density and temperature profiles 
calculated with different cuts in volume using a broken power-law 
function of eq.[T]in the logarithmic form. 

For the density profiles the internal slope does not depend 
much on the physics or the volume cutting scheme adopted: this 
shows that our method does not introduce any bias to the results of 
density for the internal regions of clusters. The external slope for 
the density profiles does not show any significant dependence on 
the physics but it decreases regularly when we include more and 
more cool and diffuse gas: this indicates that to obtain a fair repre- 
sentation of the cluster temperature in the outskirts it is necessary 
to consider the 99 per cent volume profiles. If we include all the 
gas particles the value of ^2 continues to go down, but the profiles 
begin to show features due to the cold and gas and the of the fit 
increases. For this reason and since the 99 per cent profiles show a 
remarkable regularity, as already said in the comments to Fig.|2l we 
can assume these as representative of the global cluster properties, 
although there may be some small systematic uncertainties due to 
our method. 

As already discussed in Section [3] the presence of cooling in 
the simulations creates steeper temperature profiles at the centre, 
and this can be seen by the different values of b\ in Table IA2I 
Again, like for the temperature profiles, the value of the internal 
slope is not influenced by the volume selection scheme. On the 
contrary the values of 62 are much less affected by the physics but 
show an increasing trend when adding colder gas, considering also 
a significant scatter between the 9 clusters. The difference between 
the 99 and 100 per cent values for ^2 is particularly evident with 
the ovisc model, but, again, given the higher values of the and 
the regularity of the 99 per cent profiles, we can draw the same 
conclusions as for the density profiles. 
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Table Al. Best-fit results for the inner slope hi , outer slope b2 and break radius of the gas density profiles fitted with a broken power-law. The profiles are 
computed considering different cuts in volume, from 50 to 100 per cent. 



Volume cut 


50% 


80% 


90% 95% 


99% 


100% 


ovisc 
csf 


2.47 ±0.03 
2.42 ±0.02 


2.50 ±0.03 
2.45 ±0.01 


bi 

2.50 ±0.03 2.47 ±0.03 
2.43 ±0.02 2.37 ±0.02 


2.46 ±0.03 
2.46 ±0.03 


2.55 ±0.02 
2.48 ±0.03 


ovisc 
csf 


3.86±0.55 
3.85 ±0.45 


3.78 ±0.43 
3.76±0.37 


bi 

3.71 ±0.39 3.60±0.37 
3.66±0.33 3.53 ±0.32 


3.38±0.29 
3.41±0.19 


3.12±0.32 
2.98 ±0.38 


ovisc 
csf 


1.12±0.14 
1.15±0.14 


1.20 ±0.14 
1.24 ±0.12 


RblRioo 
1.20±0.15 1.14±0.15 
1.23±0.13 1.10±0.16 


1.14±0.20 
1.31 ±0.28 


1.23 ±0.14 
1.34±0.15 



Table A2. Like Table|AT]but for the mass-weighted temperature. 



Volume cut 50% 80% 90% 95% 99% 100% 



ovisc 
csf 


0.43 ±0.03 
0.47 ±0.05 


0.43 ±0.02 
0.55 ±0.02 


bi 

0.47 ±0.02 0.48 ±0.02 
0.55 ±0.02 0.55 ±0.02 


0.41 ±0.02 
0.57 ±0.02 


0.54 ±0.03 
0.59 ±0.02 


ovisc 
csf 


1.30 ±0.63 
1.35 ±0.67 


1.22±0.53 
1.38±0.51 


bi 

1.31 ±0.54 1.49±0.51 
1.46 ±0.47 1.62 ±0.49 


1.73±0.55 
1.90 ±0.46 


2.15±0.65 
2.07 ±0.64 


ovisc 
csf 


1.21 ±0.24 
1.20 ±0.24 


1.19 ±0.25 
1.44 ±0.22 


1.33±0.21 1.40±0.17 
1.46±0.20 1.48±0.18 


1.29 ±0.24 
1.56±0.18 


1.49±0.18 
1.50±0.17 



